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Abstract. We argue that the non-detection of oscillations in Procyon by the MOST satellite reported bv lMatthews et alJ 120041) 
is fully consistent with published ground-based velocity observations of this star. We also examine the claims that the MOST 
observations represent the best photometric precision so far reported in the literature by about an order of magnitude and are the 
most sensitive data set for asteroseismology available for any star other than the Sun. These statements are not correct, with the 
most notable exceptions being observations of oscillations in a Cen A that are far superior. We further disagree that the hump 
of excess power seen repeatedly from velocity observations of Procyon can be explained as an artefact caused by gaps in the 
data. The MOST observations failed to reveal oscillations clearly because their noise level is too high, possibly from scattered 
Earthlight in the instrument. We did find an excess of strong peaks in the MOST amplitude spectrum that is inconsistent with a 
simple noise source such as granulation, and may perhaps indicate oscillations at roughly the expected level. 
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1. Introduction 

The MOST sa tellite (Microvariability and Oscillations of Stars; 
IWalker etalJl2003l) comprises a 15-cm telescope and CCD de- 
tector design ed to perf orm high-precision photometry on bright 
stars. Matthe ws et alJ (|2004, hereafter M04) reported observa- 
tions of the F5 star Procyon A made almost continuously over 
32 days. The Fourier amplitude spectrum showed no evidence 
for oscillations and M04 concluded that if there are p-modes 
in Procyon, they must have lifetimes less than 2-3 days and/or 
peak amplitudes <15 parts per million. We agree completely 
with this conclusion but we do not agree with some of the 
other statements in that paper, as outlined in the following sec- 
tions. For other discussion of the issues raised in M04, see also 
lElsworth & Thompson! ( 120041) . 



2. Lowest noise ever? 

We have measured the noise level in the MOST amplitude spec- 
trum from the data plotted in Figure 3 of M04, which was pro- 
vided to us in digital form by the MOST team. We find that the 
noise is about 2 ppm (parts per million) at frequencies above 
2 mHz and 2.5 to 4.0 ppm between 2.0 and 0.5 mHz, in agree- 



ment with the analysis by the MOST team (J. Matthews, priv. 
comm.). 

M04 described the MOST result as "the best photomet- 
ric precision so far reported in the literat ure JGilliland et alJ 
Il993t iKurtz et aljfl997t iRetter et alJEool by about an order 
of m agnitude." 1 This statement is not correct. The work by 
iGilliland et alJ Jl993l) was a ground-based photometric cam- 
paign to monitor several stars in the open cluster M67, most 
of which had spectral types similar to Procyon. Table 11 of 
that paper gives the final noise levels in the amplitude spec- 
tra in the frequency range 2.5^1.5 mHz (column 4). Values for 
the best ten stars (excluding Star 12, which is a red giant) are 
6.4-9.7 ppm. These noise levels, obtained after one week, are 
certainly higher than the 2 ppm reported for the 32-day MOST 
observations, but only by factors of 3-5 and not an order of 
magnitude. We also note that the noise level for a given ob- 
serving time is very similar for the two data sets. 



1 Due to an error in ty pesetting, the reference numbers in the text 
of Matthe w's et alJ 1200 4") do not match those in the list of references. 
The correct numbers are given in the Erratum, which we repeat here 
for convenience: References 1 to 26 should be, respectively: 1, 10-16, 
2, 3, 17-26 and 4-9. 
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Th e observations by iKurtz et all l l 19971) and iRetter et"afl 
(2004) cited by M04 targetted Ap stars and were indeed much 
less precise than the MOST data. However, the IRetter et"a"fl 
observations were made from space using the star tracker on 
the WIRE satellite, and the same instrument achieved far bet- 
ter precision (and g ood evidence for oscillations) on the solar- 
type star a Cen A JSchou & Buzasil 200 ll) . Indeed, the noise 
level from those observations was substantially below 1 ppm 
and therefore considerably better than that from the MOST ob- 
servations of Procyon. 

M04 also stated that the duty cycle of 99% for 32 days 
makes the MOST data the most sensitive data set for astero- 
seismology available for any star other than the Sun. While 
the duty cycle is certainly tremendously good, one must also 
take the precision into account and such a grand claim is 
certainly not justified. The noise at high frequencies (about 
2 ppm above 2 mHz) is about half the peak solar signal. Many 
ground-based velocity campaigns over more than a decade 
have achieved equivalent noise levels wel l below this on 
Procyon (e.g., 8cms _1 by Marti c et aDl2004l which is a third 
of the peak solar signal) and other stars (for recent r eviews, see 
iBouchv & Carrierl2003UBedding & Kieldserl2003T) . 

Recent ground-based velocity observations of the a Cen 
binary system are particularly not eworthy. Noise levels hav e 
reached 4.3 cm s -1 i n a Cen A ( Bouchv & Carrier 2002), 
3.8 cm s _I in a Cen B JCarrier & Bourb am27)03l)a nd, most re- 
cently, 2.0 cms 1 in a Cen A jButler et alJl2004 . The latter 
value is a mere tenth of the peak solar oscillation amplitude, 
and all three campaigns produced data sets that are much more 
sensitive for asteroseismology than the MOST Procyon data. 
This is reflected in the fact that they actually succeeded in de- 
tecting oscillations in stars whose amplitudes are substantially 
lower than Procyon's. 

3. Previous work on Procyon 

Several radial-velocity studies of Procyo n have detected an 
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1999; 


Bouchv et alJ 


2002; 


Kambe et al] [2003; 


Martic et al. 2004; Eggenberger et al. 2004). The reported peak 



amplitudes of 50-70cms _1 (not 1ms -1 , as quoted by M04) 
are l ower than predicted by theoretical models of oscilla- 
tions JChristensen-D alsgaard & Frand s"en1ll983t iHoudek et alJ 
1999, 1). Even before most o f these observations were made, 
Kieldsen & Beddini i 19951) concluded that F-type stars, 



namely Procyon and several members of M67, must oscillate 
with amplitudes less than has generally been as sumed. This 
point has been raised several times since then ( Martic et al. 
1999t ISamadi & Houdekl |200(A iBedding & Kieldsenl 120031 
Mar tic et al]2004 . The knowledge that Procyon oscillates with 
an amplitude only 2-3 times solar, rather than the theoretically 
predicted value of at least 4-5 times solar, has therefore been in 
the literature for many years. The MOST upper limit of 15 ppm 
(four times solar - see Sec. @J is consistent with this and is 
therefore hardly a surprise. 

It should also be noted that the theoretical predictions 
are based on our c learly incomplete u nderstanding of stellar 
convection. Indeed, Christensen-Dalsgaard & Frandsen ( 1983) 



noted that mixing-length theory overestimates the convective 
flux in relatively hot stars and that predictions of oscillation am- 
plitudes may have to be reduced. Interestingly, hydrodynami- 
cal simulations of convection for a variety of stellar parameters 
bv lStein e t al. (2004) also showed a substantial increase in the 
energy input from convection to the mode s with increasing ef- 
fective temperature. However. IStein et al] did not estimate the 
damping rate of the oscillations and hence could not provide 
estimates of the resulting amplitudes. 

Although several independent detections of excess power 
in Procyon have been made from velocity measure- 
ments, it has been difficult to obtain a power spec- 
trum that shows a clear p-mode structure. There is con- 
sensus that t he so-called large frequency separat i on is 
about 55 ^Hz jMosseretalJll998t iMartic et alJll999t l2004 
Egg enberger et alJ l2004. which agrees well with theo- 
retical models (Barbanetal. 1 19991: IChabover et alJ I l999»i 
| Pi Mauro & C hristensen-Dalsgaardl l200lUProvost et al] [2002: 
Kervella et al. 2004), but there is not yet agreement on the ac- 
tual frequencies. This is presumably due either to insufficient 
high-quality data (i.e., confusion from aliases in the spectral 
window) or to intrinsic properties of the star (i.e., short mode 
lifetimes and/or modes shifted from the asymptotic relation by 
avoided crossings). 



4. Could MOST have detected oscillations? 

For reference in the following discussion, we note that the 
strongest oscillation modes in the Sun have velocity ampli- 
tudes of 20-25 cm s _1 and intensity amplitudes in the blue- 
visual part of the spectrum of 4.5-5.0ppm. Of course, these 
peak heights decrease when the observing time becomes sig- 
nificantly longer than the typical lifetime of the modes. 2 Note 
that M04 quo ted mode lifetimes for th e Sun of days to weeks, 
referring to Toutain & Frohlicl] i 19921). However, the recipro- 
cal linewidths given in Table 6 of iToutain & Frohlichl dl992h 
must be divided by n to convert to mode lifetimes, as discussed 
in section 4 of that paper. The dominant modes in t he Sun have 
lifetimes of 2-4 days (see also Cha plin et al.l 1997). 

We must also keep in mind that intensity variations are a 
consequence of temperature variations and, while the two are 
proportional, the constant of proportionalit y depends on the ef- 
fective temperature of the star (see iKieldsen & Beddinglll995l 
Eq. 5). In a star with the same physical oscillation amplitude 
as the Sun (as measured in velocity) but with the effective tem- 
perature of Procyon (T e ff = 6530 K), we would observe an in- 
tensity amplitude at these wavelengths that is T^/T^ = 0.78 
times the intensity amplitude we observe in the Sun. Thus, the 
upper limit of 1 5 ppm in Procyon reported by M04 corresponds 
to a "physical" oscillation amplitude of 4 times solar, which 
translates to a velocity amplitude of just below lms -1 . A ve- 
locity amplitude of 60 cms 1 , as reported from ground-based 
Doppler measurements, corresponds to an intensity amplitude 
of 9.5 ppm, which is below the limit set by M04. 



2 We follow the usual definition that the lifetime of a damped oscil- 
lator is the time for the amplitude to decay by a factor of e. 
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Fig. 1. The amplitude spectrum of Procyon observed by MOST (grey), compared with simulations (black) that included oscilla- 
tions and granulation at the expected levels. The simulations were made for three different mode lifetimes, both with and without 
instrumental noise. 



To illustrate some of these points, Fig. ^ compares the 
MOST results with simulations. In each panel, the ampli- 
tude spectrum in grey is the spectrum published by M04 (the 
notches correspond to harmonics of the satellite orbital fre- 
quency, where they filtered to removed power from stray light 
artifacts). The amplitude spectra shown in black are simula- 
tions of the signal expected from Procyon if the p-mode oscilla- 
tions have peak amplitudes of 10 ppm, which matches the level 
inferred from ground-based velocity observations. We show 
simulations for three different values of the mode lifetim e (14, 
3 and 1 d), made using the met hod described by IStello et al] 
J2004 and lBedding etail J2004 . 

The simulations also included granulation noise at 1.5 
times the solar level (see below), which is responsible for the 
noise floor that rises towards lower frequencies in the upper 
three panels. In those panels, no photon or instrumental noise 
was included in the simulation. In the lower three panels, on the 
other hand, we added sufficient noise to bring the total power up 
to the level observed by MOST. As expected, we see that short- 
ening the mode lifetimes has the effect of lowering the peak 
heights, making the oscillations even harder to detect, as also 
found by M04 (Supplemental Material). It is apparent that if 
Procyon has oscillations with the amplitudes suggested by ve- 
locity measurements, we would not expect to see them clearly 
in the MOST spectrum unless the mode lifetimes were much 
longer than in the Sun. 



5. Has MOST detected granulation noise? 

M04 suggested that the increase in noise towards low frequen- 
cies that is seen in their amplitude spectrum could be caused 
by light variations due to stellar granulation. We would expect 
granulation to be much easier to detect in intensity than veloc- 
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Fig. 2. Smoothed power density spectra from photometry 
of Procyon (MOST), the Sun (VIRGO, green channel), and 
an F-type star in the cluster M67 (ground-based telescopes). 
Also shown are Balmer-line equivalent-width measurements of 
Procyon (EW, with 2-cr error bars). 



ity, but the level of non-white noise in the MOST data seems 
too high for this. 

We illustrate this in Fig. |2 which shows smoothed power 
density spectra (power per frequency resolution bin). Power 
density is the appropriate way to compare noise sources, 
whether the noise arises from granulation, photon statis- 
tics or instrumental effects, because it is independent of 
the length of the observa tions (see, e.g., Appendix A.l of 
iKieldsen & Bedding! Il995l) . The data for the Sun are full- 
disk intensity measurements from the green channel of the 
VIRG O three-channel sunp hotometer on the SOHO spacecraft 
jFrohlich et al Jl 1 997tlPalle et alJll999h . The figure shows that 
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the power density of the MOST noise is about 16 times higher 
than the solar granulation power (4 times in amplitude). This 
seems implausibly large, especially as Ba lmer-line equiv alent 
width measurements of Procyon by Kields en'et al.l \ 1999!) pro- 
duced an upper limit on granulation - and a tentative detection 
- of about twice solar in power ( 1 .4 in amplitude). Those points 
are also shown in Fig. [2] after conversion from equivalent- 
width to intensity (see Sec. 5 of that paper). We also see that the 
noise level in the MOST spectrum of Procyon is higher than the 
no ise level measured fr om photometry of F-type stars in M67 
bv lGilliland etaD dl993h . 3 

Given these comparisons, it seems unlikely that the noise 
in the MOST spectrum is due to granulation in Procyon. We 
suggest that the most likely origin of the large amount of non- 
white noise in the MOST data is instrumental, possibly from 
variations in the amount of scattered Earthlight reaching the 
detector (see M04). 

6. Could the ground-based results be artefacts? 

M04 suggested the possibility that apparent p-mode power 
identified by ground-based observers may be an artefact of in- 
sufficient gapped sampling of granulation variations. Is it possi- 
ble that the hump of excess power at 0.5-1 .5 mHz that has been 



seen repeatedly in the velocity data is 


an artefact? This was in 


fact suggested by Kieldsen & Bedding 


( 1995) as a likely cause 


for the power excess found by Brown et alJ dl 99 ll). on the ba- 



sis that the observations were high-pass filtered. Such a filter, 
when applied to a Fourier spectrum dominated by non-white 
noise from instrumental effects, will indeed produce a hump. 
However, subsequent repeated independent confirmations of 
the power hump, most of whic h were high-pass filtered m uch 
less strongly (or not at all), led Beddi ng & Kieldseiil d2003l) to 
acknowledge that this excess has a stellar origin. 

We also note that in discussing Fig. 7 from IMartic et alJ 
(1999), M04 have chosen a figure showing only a subset of 
the data discussed in that paper, which was deliberately cut to 
allow comparison with observatio ns of another star (n Cas). On 
the other hand, Figs. 4 and 8 from M artic et alJdl999t) were not 
high-pass filtered at all and clearly show the excess power. 

The suggestion by M04 is that non-white noise from stel- 
lar granulation could have produced a hump of power because 
the ground-based observations have daily gaps. While gaps in 
the observing window certainly introduce false peaks (aliases), 
these are always symmetrical about the original peak and are 
very close to it (small multiples of + 1 1 .6 //Hz). It is not possible 
using this mechanism to turn a monotonically varying power 
spectrum (as produced by granulation or, indeed, by most in- 
strumental effects) into a hump of excess power, as we have 
verified using simulations. We therefore disagree strongly that 
the excess power seen by ground-based observers can be ex- 
plained as an artefact introduced by gaps in the observations. 

3 Note that all four power density spectra in Fig. [2] were scaled in 
the same way, such that a sinusoid with an amplitude of 1 ppm ob- 
served for 10 6 s gives a peak in power density of height 1 ppm 2 //jHz. 
This explains the factor of t wo between the VIR GO power density in 
the figure and that plotted bv lPalle etal1 fl9991. who used a different 
scaling convention. 




Amplitude (normalised) 

Fig. 3. Distribution of peaks heights in the MOST amplitude 
spectrum of Procyon, compared with pure noise and also with 
three frequency bands in the GOLF solar spectrum (see text). 

Fin ally, we note that in response to M04, iBouchv et alJ 
(2004) have published a short sequence of very high-precision 
velocity measurements of Procyon obtained with the HARPS 
spectrograph. These show clear periodicities with typical peri- 
ods of 18 minutes and an excess of power in the Fourier spec- 
trum around 1 mHz, in perfect agreement with the previous ve- 
locity measurements. R ecent evidence for t he power excess has 
also been presented bv lClaudi et a l. (2005) on the basis of ob- 
servations with a longer timespan but lower precision. 

7. Evidence for oscillations in the MOST data? 

We have examined the MOST amplitude spectrum for evidence 
of oscillations. Based on the ground-based observations, oscil- 
lations could be present at low signal-to-noise. Indeed, inspec- 
tion of Fig. 3 in M04 shows a few peaks that are well above the 
surrounding noise and it is possible that some of these are due 
to oscillations. 

Comparison with simulations shows that there are many 
more strong peaks in the amplitude spectrum than would 
be expected from noise. To quantify this, we examined 
the cumulative distribution of peak heights, a method 
that has previously been appli ed in both helioseismology 
(e.g.. IPelache & Scherreil fl983l) and asteroseismology (e.g., 
IBrown et alJIl 99 lit . We first removed the overall slope by di- 
viding by a smoothed version of the spectrum. This flattens the 
spectrum and allows us to make meaningful counts of the num- 
ber of points above a given amplitude. The resulting distribu- 
tion is shown in Fig. [3] where the horizontal axis is amplitude 
divided by the mean amplitude. 

On the same graph we show the distribution resulting from 
noise, which we measured by averaging a large number of sim- 
ulations. It is clear that the MOST spectrum has an excess of 
strong peaks. For example, there are more than ten times as 
many peaks with height > 3.5 times the mean than expected 
from noise. 

Also shown in Fig. [5] are distributions of peak heights for 
the Sun, based on a 30-day series of full-disk velocity observa- 
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tions take n by the GOLF instrument on the SOHO spacecraft 
( Ulri ch et al.l20 0C0. We show results for three frequency bands 
in the solar amplitude spectrum. Band 1 (0.4-1. 4 mHz) is at 
frequencies below the solar oscillations and contains only noise 
from granulation. We see that this distribution closely matches 
our noise simulations. Band 2 (1.4-2.4 mHz) contains granula- 
tion noise plus the weak low-frequency solar p-modes. Those 
modes have long lifetimes compared to the length of the time 
series and so the distribution contains unresolved high peaks. 
Band 3 (3.8^4.8 mHz) contains the weak highest-frequency so- 
lar p-modes, which have short lifetimes. Thus, this band has 
many resolved modes, which produces a broader distribution 
of peaks. 

Comparing the different curves in Fig.0 we conclude that 
the distribution of peak heights in the MOST amplitude spec- 
trum is not consistent with being solely due to noise from gran- 
ulation or similar sources. Rather, there is evidence for an ex- 
cess of strong peaks that is consistent with the presence of os- 
cillations at approximately the level expected. However, the un- 
usual peak distribution could arise from other causes (as turned 
out to be the case for th e solar measurements analysed by 
iDelache & ScherreJll983l) . In particular, the filtering by M04 
to remove power at the harmonics of the satellite's orbital fre- 
quency may have distorted the peak-height distribution and fur- 
ther analysis of the original time series is needed to help decide 
this question. 

8. Conclusions 

Our main conclusions are as follows: 

1. The upper limit placed on oscillations in Procyon by the 
MOST satellite is consistent with previous ground-based 
velocity measurements and with the fact that Procyon os- 
cillates with an amplitude significantly below theoretical 
predictions, as has been discussed in the literature for many 
years. The conclusions by M04 that their results "defy ex- 
pectations from the Sun's oscillations and previous the- 
oretical predictions" and that "target selection for future 
planned asteroseismology space missions may need to be 
reconsidered, as will the theory of stellar oscillations" are 
overstated. 

2. The claims by M04 that the MOST data represent the best 
photometric precision so far reported in the literature by 
about an order of magnitude and are the most sensitive data 
set for asteroseismology available for any star other than 
the Sun are both inaccurate. For example, published obser- 
vations of oscillations in a Cen A in both photometry and 
velocity are far superior. 

3. We find it unlikely that the non- white noise in the MOST 
amplitude spectrum is due to stellar granulation. Instead, 
the most likely explanation lies in variable amounts of scat- 
tered Earthlight reaching the detector. This instrumental 
noise is probably the reason that oscillations could not be 
detected in Procyon. 

4. The distribution of peak heights in the MOST amplitude 
spectrum is not consistent with a simple noise source such 
as granulation. The excess of strong peaks may indicate os- 



cillations, but further work is needed before drawing a firm 
conclusion. 

The excellent observing duty cycle of MOST makes this a 
potentially powerful instrument for asteroseismology. It is pos- 
sible that more sophisticated treatment of the stray light will 
reduce the noise level enough to reveal a clear detection of 
oscillations in Procyon. If the non-white noise cannot be re- 
duced, the best targets will be those with larger amplitudes than 
Procyon, such as more evolved stars (e.g., rj Boo) and classi- 
cal pulsators such as 6 Scuti stars and roAp stars. MOST has 
tremendous potential to contribute greatly to the study of these 
types of stars. 
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